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Abstract 

Recently, as demonstrated by an antiferromagnetic spin-lattice application, 
we have successfully extended the coupled-cluster method (CCM) to a varia- 
tional formalism in which two sets of distribution functions are introduced to 
evaluate Hamiltonian expectation. We calculated these distribution functions 
by employing an algebraic scheme. Here we present an alternative calculation 
based on a diagrammatic technique. Similar to the method of correlated-basis 
functionals (CBF), a generating functional is introduced and calculated by a 
linked-cluster expansion in terms of diagrams which are categorized and con- 
structed according to a few simple rules and using correlation coefficients and 
Pauli exclusion principle (or Pauli line) as basic elements. Infinite resumma- 
tions of diagrams can then be done in a straightforward manner. One such 
resummation, which includes all so-called ring diagrams and ignores Pauli 
exclusion principle, reproduces spin- wave theory (SWT). Approximations be- 
yond SWT are also given. Interestingly, one such approximation including 
all so-called super-ring diagrams by a resummation of infinite Pauli lines in 
additional to resummations of ring diagrams produces a convergent, precise 
number for the order-parameter of the one-dimensional isotropic model, con- 
trast to the well-known divergence of SWT. We also discuss the direct relation 
between our variational CCM and CBF and discuss a possible unification of 
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the two theories. 
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I. INTRODUCTION 



A microscopic quantum many-body theory is mainly to study correlations between the 
constituent particles of a quantum system. One of the most successful quantum many-body 
theories is the method of correlated basis functionals (CBF) [1] which includes real-space 
many-body correlation functions in the ground state and employs similar techniques as in 
classical statistical mechanics to calculate the corresponding distribution functions. The 
CBF has proved to be one of very few many-body theories capable of dealing with strongly 
correlated boson systems in liquid phase. Another successful quantum many-body theory 
is the coupled-cluster method (CCM) [2-4] in which excitation operators with respect to 
an uncorrelated model state are employed to construct the many-body correlations in an 
exponentiated operator factor in the ket ground-state. The distribution functions are not 
needed in the CCM as the Hamiltonian expectation value is straightforwardly calculated as 
a finite order polynomial in terms of the correlation coefficients. This is because the bra 
ground-state in the CCM is not the hermitian conjugate of the ket ground-state but is in a 
simple, linear form [5]. The CCM has proved to be one of most powerful techniques in cal- 
culating ground-state energy for many non-liquid fermion systems such as atoms, molecules 
and electron gas [6]. However, a general many-body theory capable of dealing with strongly 
correlated fermion systems in liquid phase is still needed. 

In our earlier papers [7], by an application to bipartite quantum antiferromagnetic lattice 
systems, we have extended the CCM to a variational formalism where, in contrast to the 
traditional CCM, the bra and ket ground-states are hermitian conjugate to one-another. 
Two sets of the distribution functions were introduced to evaluate Hamiltonian expecta- 
tion. We developed an algebraic scheme to calculate these distribution functions through 
self-consistent sets of equations. In this paper, we present an alternative scheme based on 
diagrams to calculate these distribution functions. As in the CBF, a generating functional 
is introduced and calculated by a linked-cluster expansion in terms of diagrams. These dia- 
grams are constructed according to a few simple rules and using only three basic elements: 
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(a) dots representing the ket-state correlation coefficients, (b) exchange lines representing 
the bra-state correlation coefficients, and (c) Pauli line representing Pauli exclusion princi- 
ple manifested by the spin-1/2 operator property (s^O 2 = in our spin model example. In 
this fashion, infinite diagrams can be resummed in a straightforward manner. As a simple 
application in our spin model, the spin-wave theory (SWT) [8,9] is reproduced by including 
all so-called ring diagrams without any Pauli line. Approximations beyond SWT by resum- 
mations including Pauli lines are also given. One such approximation, which includes all 
so-called super-ring diagrams by a resummation of infinite Pauli lines in addition to resum- 
mations of ring diagrams, produces a convergent, precise result for the order-parameter of 
the one-dimensional isotropic model, contrast to the well-known divergence of SWT. Fur- 
thermore, the diagrammatic analysis discussed here forms a basis for a possible combination 
of the variational CCM and CBF. Such a unified theory may prove to be capable of dealing 
with strongly-correlated fermion system in liquid phase. 

II. THE VARIATIONAL COUPLED-CLUSTER METHOD 

We take the spin-1/2 antiferromagnetic XXZ model on a bipartite lattice as our example. 
The Hamiltonian is given by 

H = \ E H l,l+n = \Y. ( As l S l+n + ^V+n + \ S T 4+n) > W 
l,p l,n 

where A > is the anisotropy constant, the index / runs over all lattice sites, n runs 
over all z nearest-neighbour sites, and s ± are the usual spin raising (+) and lowering (— ) 
operators. As in our earlier work using the traditional CCM [10], we take the model state 
as the classical Neel state with alternating spin-up and spin-down sublattices. As before, 
we shall exclusively use index % for the spin-up sublattice and the index j for the spin-down 
sublattice. The many-spin correlations in its ground state of Eq. (1) can then be included by 
considering the excited states with respect to the uncorrelated model state. These excited 
states are constructed by applying the so-called configuration creation operators C\ to the 
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Neel model state with the nominal index / labelling these operators. In our spin model, the 
operators C\ are given by any combination of the spin-flip operators to the Neel state s~ 
and s~j~; the index / in this case corresponds to the collection of the lattice indices (i's and 
j's). The hermitian conjugate operators of Cj are the configuration destruction operator 
Cj, given by any combination of s+ and sj. For example, the two-spin flip creation operator 
is given by Cjj = s~Sj~/2s, and their destruction counterpart, Cy = sfs~/2s, where s is 
the spin quantum number. Although we are mainly interested in s = 1/2 in this article, we 
keep the factor of l/2s for the purpose of comparison with the large- s expansion. 

As discussed in details in our earlier paper [7], we use Coester representation for both 
the ket and bra ground-states and write 

|*> = e 5 |$>, S = J2FjCj; (¥| = ($ |e* S^FA (2) 
i i 

with 

N/2 - - + + 

EW) = E E h,-,n... n J \ (3) 

for the ket state and the corresponding hermitian conjugate of Eq. (3) for the bra state, 
using notation Fj = fi 1 ...^ 1 ... for the bra-state coefficients. The coefficients {Fj, Fj} are then 
determined by the usual variational equations as 

5Fj SFj ' w (^|^) v ; 

We define the so-called bare distribution functions as 

9i = (Ci), 9i = (Cj), (5) 

where we have exchanged the definition of gj with that of gj as compared with those in 
Ref. 7 for purely notational reason. The Hamiltonian expectation (H) is shown, in general, 
to be a function containing up to linear terms in gj and gj and finite order polynomial in 
Fj (or in Fj) in Eq. (21) of Ref. 7: 

(H)=H(g I ,g I ,F I )=H(~g I ,g I ,F I ). (6) 



Two systematic schemes have been developed for calculating the distribution functions of 
Eqs. (5): one is algeraic and the other is diagrammatic. In the algebraic approach, by taking 
the advantage of the properties of the operators, it is straightforward to derive the following 
self-consistent sets of equations for the distribution functions 

gi = G(gj,Fj), ~ 9l = G( 9j , Fj) , (7) 

where G is a function containing up to linear terms in gj (or gj) and finite order polynomial 
in Fj (or Fj). Eqs. (7) are solved for gi and gj as a function of Fj and Fj. The variational 
Eqs. (4) are then carried to determined the optimum Fj and Fj. In this algebraic calculation, 
direct comparison with the traditional CCM can be made. It is shown that the CCM 
is a linear approximation to one set of distributions as simply g I m Fj, which is a poor 
approximation for the spin-spin correlation function and low-lying excitations. More detailed 
comparison was given in Ref. 7. Here we present the diagrammatic scheme similar to that 
in CBF to calculate these distribution functions. We like to point out that Eqs. (2) and 
(4)-(7) are the main general equations of the variational CCM. 

III. DIAGRAMMATIC REPRESENTATION OF GENERATING FUNCTIONAL 

In this section, we calculate the bare distribution functions gj and gj of Eq. (5) by 
employing a diagrammatic scheme. As a demonstration, we consider a simple truncation 
approximation in which the correlation operators S and S of Eqs. (2) retain only the two-spin 
flip operators as (the so-called SUB2 approximation as defined in Ref. 10), 

ij ij ij ij 

Using the usual angular momentum commutations [sf,s^] = ±sf5«', [s^,s ; 7] = 2sf5jj/, and 
the Neel state eigenequations, sf|3>) = s|$), s||$) = — s|<&), it is a straightforward calcula- 
tion to derive expectation value of any physical operators in terms distribution functions of 
Eqs. (5). In this approximation, for example, the expectation value of Eq. (1) is given by 
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(ify) =A{s z i s z j ) + -(g ij + ~g ij ), (9) 



where (s*Sj) is calculated as 



( sZ i s j) = -s 2 + s \J2pi'j + J2pif - (J2pij',i'j + pij) , 

Pij is the usual full one-body distribution function defined as 



(10) 



Pij = fijdij = fij ^ j (11) 

and where Pij^y is the full two-body distribution function define as 

Pij,i'j' = fijfi'j'9ij,i'j' = fijfi'j' (2s) 2 ' 

The order parameter is given by 

(st) = s-p, (13) 

where p = J2j Pij, taking the advantage of translational invariance. 
We define a generating functional W in the usual fashion as, 

W = \n(^), (14) 

so that the bare and full distribution functions can be simply expressed as functional deriva- 
tives of W. For example, the one-body and two-body bare functions are given by 

9i = (Cl) = §- , 9u = (C\CI) = P^ + hh , (15) 

where, for simplicity, we have employed notation 1 = so that f\ = fi 1 j 1 etc.; and the 

structure function S 12 has the usual relation as in the CBF as 

S12 = fiTT = Pi^i2 + P12 - P1P2 , (16) 

where p\ = Pi 1 j 11 etc. 

We now write W in terms of a linked-cluster expansion 



W = sum of all linked cluster contributions . (17) 

The main task of this section is to find a diagrammatic scheme to categorize this expansion. 

We first expand the ket-state operator in the simplified notation, e s = 1 + S + ^S 2 H = 

1 + f\C\ + ^fif 2 ClC 2 + ' ' '5 where in the last equation, the summation over all indices is 
understood. The normalization integral, 

= 1 + f v h{C v C\) + ^f 2l f vfl f 2 (C 2l C v Clcl) + • • • , (18) 

can be evaluated straightforwardly for the first few terms. In the above series, the primed 
indices are used for bra state expansion. We notice that each term of Eq. (18) contains equal 
number of creation and destruction operators (otherwise, the expectation is zero). 

The first-order expectation is easily calculated as (CyC\) = (2^($l s 7 1 , s ^Av s / 1 1^) = 
Si 1 ,i 1 Sj i ,j 1 . Hence we have, writing out the summation explicitly, 

1st order = £/i/i . (19) 
1 

The calculation of the second-order expectation (C 2 'Ci'C\C 2 ) is slightly more complicated. 
We first consider the case of 1 7^ 2 (i.e., %\ 7^ % 2 and ji 7^ j 2 ). There are four nonzero terms 

(^i^il^/h + ^i 1 ,i2^i 2 'il) (j>3i'3l ^3 2 '32 ^3i'32^j 2 rjl) 
= ( ( ^i 1 /ii^i 2 /i2 + ^i 1 ,i2^i 2 'il) {^3i'3l^3 2 ih + ^3i'32^j 2 rji) ■ 

The cases when i x = i 2 and/or j\ = j 2 can be easily accounted for by introducing a factor 
involving the usual delta functions as, 



with a definition, 



A i2 = -tt^t (S ilia + S jlja ) + — — 2 5 ili2 5 jlh . (20) 
(2s) (2s) 

This is because (sj 7 ) = (ft^) = ^ or s = a manifestation of Pauli exclusion principle. 
The second-order contribution is hence derived as 
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(2!) 



1 r 

2! - 



where the second term inside the square brackets clearly represents the so-called exchange 
contributions. We notice that we did not consider explicitly the Pauli exclusion principle 
for the bra-state operators in the above derivations as the delta functions for the ket state 
operators also take this principle into account due to the fact that each of the bra-state 
operators always need to match one of the ket-state operators in order to give nonzero 
contribution. This is also true for higher-order terms. We also notice the expression of 
Eq. (20) is in fact also correct for spin quantum number s > 1/2 because, for a general s, 



etc. For higher-order terms in the expansion of Eq. (18), the extension of Pauli exclusion 
principle can be simply written as a product of two-body factors as 



We notice that the above product in general is not exact any more but an approximation 
for s > 1/2 as the three-body effects (e.g., from (s^) 3 when i\ — i 2 — 23) have been ignored. 

In order to extend to higher-order calculations including the exchange contributions, we 
need a systematic graph representation. For this purpose, as shown in Fig. 1, we use a solid 
dot to represent the ket state coefficient fi with 1 = (iiji) as defined earlier; a (directed) 
exchange line drawing from i\ to j 2 to represent the bra state coefficient fi^; and a Pauli 
(dashed) line drawing between any two dots to represent delta function A i2 of Eq. (20). 
With these graphic notations, a linked contribution is represented by a connected diagram. 
After the detailed calculations up to 5th order, we have established the following simple and 
complete rules for construction of these diagrams in the normalization integral of Eq. (18): 

• The /cth-order contribution consists of all possible diagrams involving k dots; 




n (i + A nm ) . 



(21) 



n>m 
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• In each diagram the number of dots equal to number of exchange lines; 

• A dot is always connected by exchange lines (leaving and coming) hence exchange lines 
always form loops; 

• Between any pair of dots one can draw at most one Pauli line; 

• The contribution of each diagram is divided by its symmetry factor; 

• Summations over all indices involved. 

We first consider the case without any Pauli line, A nm = 0. (This is equivalent to 
turning spin operators to boson operators as will be shown later.) For example, the first- 
order contribution of Eq. (19) is simply a dot with an exchange line leaving and coming as 
shown as diagram a in Fig. 2, where the direction of exchange line is clockwise as in most 
other diagrams (we therefore do not show arrows of exchange lines explicitly). The second- 
order contribution with A 12 = is given by two diagrams b and c in Fig. 2, namely ^(b + c). 
The 3rd-order contribution is calculated as |, (a + 36 + 2c) and is shown in Fig. 3, where 
factor 3 for diagram b is due to the three equivalent diagrams by rotation and the factor 2 
for diagram c comes from the two equivalent diagram with opposite directions, one clockwise 
the other counter-clockwise (this is referred to as parity symmetry). In similar fashion one 
can write down the 4th-order contribution as shown in Fig. 4 for the corresponding diagrams 
as 



where the coefficient numbers are the symmetry factors of the corresponding diagrams. For 
example, the factor 6 for diagram e is due to the fact that there are three equivalent diagrams 
each with parity symmetry factor of 2. The 5th-order contributions include 7 independent 
diagrams, as shown in Fig. 5, namely 




(22) 



5th order = — (a + 106 + 20c + 15d + 30e + 20/ + 24g) . 



(23) 
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We notice that, in all these results, the last term represents a ring diagram with k dots in 
the fcth order contribution. We use to represent this ring diagram with the symmetry 
factor (k — l)\/k\ = 1/k. For example, the 4th-order ring contribution is, writing out the 
summations explicitly 

^4 = 7 X! hfiihhfimhfizjtftfuh ■ (24) 

4 1,2,3,4 

Furthermore, the other terms in these fc-order contributions are simply a product of smaller 
ring contributions. This property can be extended to higher order. (For this purpose one 
needs to apply symmetric group S n to count the number of diagrams. See, for example, 
Ref. 11). We are now in position to write all contributions without any Pauli line in terms 
of these ring diagrams. The normalization integral of Eq. (18) is then written as 

<*!»>*- « = t~( R ^wr---^r 

= exp(i?i + R 2 + R 3 + • • •) • 
The corresponding generating functional W without any Pauli line is simply 

oo 

W' = W _ =J2 R k- ( 25 ) 

A nm — k=l 

To include Pauli lines (i.e. A nm ^ 0), we use notation Lk to represents the contribution 
of all linked /c-clusters and write 

W = ln(y\y) = L x + L 2 + L 3 + ■■■ . (26) 

Using the simple rules discussed earlier, without much difficulty, we can list all /c-cluster 
contributions of Lk in terms of a ring diagram Rk plus all possible ways of drawing Pauli 
lines between any pair of k dots of rings, including those pairs of dots between rings and 
those pairs of dots inside rings. In Fig. 6 we list all 3rd-order contributions in L 3 except 
Ri, Ri and R 3 . 
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IV. DIAGRAM RESUMMATIONS, SPIN- WAVE THEORY AND BEYOND 



We first consider all diagrams without any Pauli line, namely all the ring diagram con- 
tributions Rk with k — 1, 2, • • •, and show that the spin-wave theory is thus reproduced. As 
can be seen from Eq. (20), these ring diagrams represent the first order approximation in the 
large- s limit. In fact, in this limit, operators s~ and s+ behave like bosons as s~ — > \/2s a\, 
s+ — > \/2s 6] [8,9]. The corresponding wavefunction by Eq. (8) becomes the spin- wave 



3 v 3 

function as 



|*> - \* sw ) = exp E fiMb] ] |$) = J] exp (f q a\bl q ) |$> , (27) 



where the Neel state |<3>) should be considered as the vacuum state for the two sets of bosons 
a\ and &j and where, in the last equation, we have made Fourier transformations using the 
translational symmetry as, 



2 ~ 

f. . — \ " p -*q-(rj— ri) /■ 

q 

with summation over q restricted to the magnetic zone. The normalization integral of 
Eq. (27) can be easily calculated as the wavefunction is uncoupled in g-space. Using expan- 
sion exp(f q a q b l '_ q ) = En(/g a ^-?) n / ' n ' an d a simple algebra ($|a™(aj) n |$) = n!, we have the 
following well-known result (see, for example, Chapter 2 of Ref. 12), 

0M^> = II— Vt • (28) 

The corresponding generating functional is hence given by, 



w sw = H* sw \v sw ) = - Ein (i - /,/,) = E + + ■■■ 

q q L 

which is precisely the result of Eq. (25) after Fourier transformation, namely 



(29) 



W = W sw . (30) 
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Distribution functions without any Pauli line can be easily calculated using functional 
derivatives of Eq. (15)-(16) with diagrammatic representation. For example, the one-body 
bare distribution function, g[ = SW'/Sfi, is simply represented by Fig. 7, where the action of 
partial derivative is equivalent to unfolding the ring. Writing out the summations explicitly, 
we have the expansion of Fig. 7 as, 

9i — fi + X fhhhfhji + X fiihhfithhfhh ^ > (31) 

2 2,3 

and similar expansion for g[. A close inspection of g[ and g[ expansions yields self-consistency 
equations as 

9i — h + X fiihfhfiih i 9i = fi + X fiihti'ifizh > (32) 

2 2 

agreed exactly with Eq. (31) of Ref. 7 in this SWT approximation. The two-body functions 
in this approximation can also be easily obtained in this fashion as given in Ref. 7. The 
spontaneous magnetization of Eq. (13) is given by (sf) = s — p', with p' given by 



where we have used the reproduced SWT results of Ref. 7, 



(33) 



f q = f q = -Wl-H/A 2 -l), 7 9 = -E^ rn , (34) 

7g z n 

where z is the coordination number and n is the nearest-neighbour index of the bipartite 
lattice. For one-dimensional (ID) model at isotropic point A — 1, the integral of Eq. (33) 
diverges, contrast to the well-known exact result of p = 1/2 for s = 1/2 by Bethe ansatz 
(see Ref. 10 for references). 

To go beyond SWT, we need to include Pauli lines. Using the similar resummation tech- 
nique as discussed above, we express the expansion of bare one-body distribution function 
<7i in terms of diagrams as shown in Fig. 8, similar to the expansion in Chap. 9 of Ref. 12 
and in Ref. 13, after multiplying fi on both sides of the equation, 

5W 

Pi = fih — /lTT — Fig. 8 , (35) 
ofi 
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where we have done all resummations of ring diagrams as in Eq. (31) and hence all exchange 
lines in the diagrams of Fig. 8 are now function g'-, not the original exchange line function 
fij. For a simple approximation, we consider the first two diagrams of Fig. 8 as 

(36) 

2 

We notice that, after ignoring the higher-order term, Eq. (36) (without the common fi 
factor) agrees with the expression (sfsj) = (s~s^) of SWT in Ref. 9 and with Eq. (31) of 
our earlier paper Ref. 7 (after changing the sign of both 5th and 6th terms in the equation as 
they were typos). After summing over index ji with p = Phjn we have, using Eq. (20) 
for A 12 , 

" = "'-|W 2+ (2^SH) 2 . (37) 

For s = 1/2 and isotropic point A = 1, we obtain p m 0.127 for the square lattice and 0.067 
for the cubic lattice. They are smaller than p' = 0.197 and 0.078 of SWT respectively. This 
is not surprising because SWT is known to have over estimated the quantum fluctuations. 
The best numerical values for the square lattice vary from p = 0.16 to 0.19, including results 
from extrapolation of high-order localized CCM calculations [14]. For the ID isotropic 
model, however, p of Eq. (37) diverges as p' diverges as mentioned earlier. 

We next consider an approximation involving higher-order Pauli lines by including all 
higher-order diagrams similar to that of Eq. (36), as shown in Fig. 9. This infinite series can 
again be resummed in a closed form as a self-consistency equation, equivalent to replacing 
p' 2 in Eq. (36) by p2 itself as 

Pi=Pl+P / iE A i2P2. (38) 

2 

The resummation in Eq. (38) is similar to the resummation of rings in Eqs. (31)-(32), we 
therefore refer it as super-ring resummation. The numerical results for p thus obtained at 
the isotropic point for high dimensions improve slightly, as p = 0.145 for the square lattice 
and 0.068 for the cubic lattice. However, for the ID isotropic model, Eq. (38) produces a 
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convergent, precise number p = 1/2. This is interesting indeed, as the divergence of SWT has 
troubled theorists for many years. It is worth mentioning that the traditional CCM SUB2 
approximation [15] also produced a convergent result for the ID model but at A = 0.373, 
not at the isotropic point A — 1. We leave more discussion to the following section, and 
leave detailed calculations including other higher-order terms and resummations in two-body 
function p i2 and structure function Su of Eq. (16) somewhere else. 

V. DISCUSSION 

In this article, we present a diagrammatic scheme for the calculations of distribution 
functions of the variational CCM, as an alternative to the algebraic scheme published in our 
earlier papers [7] . The results of SWT are reproduced by an approximation which resums all 
ring diagrams without any Pauli line. Approximations beyond SWT can also easily be made 
by including diagrams with Pauli lines. One such approximation, which includes all super- 
ring diagrams by a resummation of infinite Pauli lines in addition to resummations of all ring 
diagrams, produces a convergent, precise number for the order parameter of the ID isotropic 
model, contrast to the divergence of SWT. This cure of SWT divergence is also interesting to 
2D models (including square and triangle lattices) as naive higher-order calculations within 
the framework of SWT are also likely to produce divergent results, despite the fact that 
the first order results are reasonable. We believe that similar resummations of super-ring 
diagrams as Fig. 9 and Eq. (38) may provide a solution for such divergent problems. We 
leave more detailed calculations to somewhere else. 

It is also possible to include in the ground state higher-order many-body correlations 
such as 4-spin-flip operators, in additional to the 2-spin-flip operators of Eqs. (8). Fur- 
thermore, as demonstrated here by the diagrammatic approach, a direct link between our 
variational CCM and the powerful CBF has now been established, as both rely on deter- 
mination of distribution functions through functional derivatives of a generating functional. 
In particular, as given by Eq. (13), particle density p in CBF is equivalent to the order 
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parameter of our spin models as (sf ) = s — p. Its diagrammatic expansions in two theories 
are similar (see Chap. 9 of Ref. 12 and Ref. 13 for more CBF details). For 2D and 3D 
lattice models, the values of density p are small compared with s. Such spin systems can 
therefore be described as dilute gases (dilute gases of quasiparticle magnons of spin waves). 
For the isotropic ID model, density p is saturated, corresponding to the order parameter 
equal to zero, a critical value. Our approximation including a resummation of super-ring 
diagrams is capable of reproducing precisely such number. It is also interesting to know that 
our diagrammatic analysis of the variational CCM is for the translational invariance lattice 
system while similar analysis in CBF is for inhomogeneous systems [12,13]. 

Clearly, for more accurate results in general, we need to include correlations between 
those quasiparticles in our ground state, and the CBF is well known to be one of most 
effective theories for dealing with such particle correlations (even when they are very strong 
as in a Helium-4 quantum liquid [1]) by systematic calculations of the important two-body 
distribution functions. We therefore propose a unified trial wavefunction \^u) as > including 
a generalized Jastrow correlation operator S° involving quasiparticle density operator s z , 



where {/,°} are the new additional variational parameters and \^f) is our variational CCM 
state of Eq. (2). The diagrammatic scheme as discussed in this article is useful for calculating 
the expansion of the new generating functional of Eq. (39). We have made progress in such 
calculations and wish to report results soon. We also believe such a unified many-body theory 
may prove to be capable of dealing with strongly correlated fermion systems in general. 
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FIGURES 



Fig. 1. Three basic elements for construction of diagrams, where simplified index nota- 
tions 1 = etc. are used. 

Fig. 2. First- and Second-order contributions in the expansion of Eq. (18). 
Fig. 3. Diagrams of third-order contributions in the expansion of Eq. (18). 
Fig. 4. Similar to Fig. 3 but for the 4th-order contributions. 
Fig. 5. Similar to Fig. 3 but for the 5th-order contributions. 

Fig. 6. Diagrams of up to third-order contributions to the generating functional W of 
Eq. (26) except ring diagrams R±, R 2 , -R3. The corresponding symmetry factors are, in the 
same order as the list of diagrams, (1/2, 1/2, 1, 1, 1/2, 1, 1/2, 1, 1/3!, 1/2, 1/3). 

Fig. 7. The ring expansion of the one-body bare distribution function gi of Eq. (31). 

Fig. 8. First few contributions to the full one-body distribution function p\ of Eq. (35), 
where open dots indicating no summations over its indices while solid dots indicating such 
summations as before. 

Fig. 9. Super-ring diagram expansion, similar to the ring diagram expansion of Fig. 7 
but now involving Pauli lines and with resummations of ring diagrams already carried out 
in all exchange lines. See texts for more details. 
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